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Linear conductance through a quantum dot is calculated under a finite magnetic field using 
the modified perturbation theory. The method is based on the second-order perturbation theory 
with respect to the Coulomb repulsion, but the self-energy is modified to reproduce the correct 
atomic limit and to fulfill the Friedel sum rule exactly. Although this method is applicable only 
to zero temperature in a strict sense, it is approximately extended to finite temperatures. It 
is found that the conductance near electron-hole symmetry is suppressed by the application of 
the magnetic field at low temperatures. Positive magnetoconductance is observed in the case of 
large electron-hole asymmetry. 
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To date, numerous theoretical as well as experimental 
studies have been carried out on the tunneling trans- 
port of an_electron through a quantum dot connected to 
the leads!? BP The combination of single electron charg- 
ing and energy level quantization causes the Coulomb 
blockade or Coulomb oscillation phenomena, which re- 
flects the particle character of an electron. Further- 
more, a dot with an odd number of electrons resembles a 
magnetic impurity coupled to the conduction electrons 
in a metal, and hence, the Kondo-type phenomenon, 
which reflects the wave characterxif electrons, has been 
predicted by several researchersBl3# and several obser- 
vations, oi. Kondo-assisted tunneling were reported re- 
cent lyBOBtP The Kondo effect in a quantum dot system 
has brought up new and interesting issues for physics, 
e.g., the tunable Kondo effect or the nonequilibrium Kondo 
effect. 

Several theoretical methods have been devised to ex- 
plain Kondo-type transport through a quantum dot us- 
ing the impurity AndeEsoa. model: second-order pertur- 
bationJiicDry (SOPT) JBHH modified perturbation the- 
ory,^ otil 1 slajifi hosnn methods noncrossing approxi- 
mation (NCA),BBE3) quantum Monte Carlo ZQMC)EP 
and rnopierical renormalization group (NRG). EH Yeyati 
et alJ& proposed an interpolative scheme that repro- 
duces the correct atomic limit in addition to the weak 
correlation limit, and applied it to transport through a 
dot. We have reinvestigated their method and modified 
it to create a more natural scheme by introducing the 
effective energy level within a dot to fulfill the Friedel 
sum rule exactlyO* The main features of a dot were suc- 
cessfully calculated by this method. Our scheme, how- 
ever, was applicable only to zero temperature in the strict 
sense of the Friedel sum formula. 

In this letter, we approximately extend our previous 
study to finite temperatures, and show the temperature 



and external magnetic field dependence of the linear con- 
ductance through a quantum dot in the Kondo regime. 
We also investigate the spin susceptibility and specific 
heat of the dot, and compare them with the Bethe Ansatz 
solution of the Anderson model to confirm the accuracy 
of our method. As a result of numerical calculation, 
the linear conductance through a dot decreases largely 
with increasing temperature near electron-hole symme- 
try. This result is consistent with those of previous theo- 
retical calc^lationsoLil'O and agrees with experimental 
results.D'EfLF As regards the magnetic field dependence, 
wc demonstrate the reduction of the conductance near 
electron-hole symmetry which is similar to the tempera- 
ture dependence. Moreover, our calculation predicts the 
positive magnetoconductance in the electron-hole asym- 
metric case. Although it was already pointed out in our 
previous paper for a limited casejls the behavior of the 
conductance is clarified in the present letter for a wider 
range of parameters. 

The Hamiltonian for a quantum dot connected to the 
leads is written as 

H = ^ Zk n k,a + e ° n<T + Un T n l 

V.k,(T (7 

where £o and U represent an energy level and the Coulomb 
repulsion in a dot. The change of Eq is equivalent to the 
change of the gate voltage in measurements. e v k denotes 
the conduction electron energy in the lead v(=R and 
L). V k v denotes the coupling between leads and a dot. 
We neglect orbital degeneracy and the k dependence of 
V k . We assume the quantum limit and strong coupling, 
therefore, the relation Ae ^> U > T ^> T should be sat- 
isfied, where Ae is the spacing of the discrete level, T is 
the resonant level width and T is the temperature. 
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Green's function for the electron in the dot is given by 

1 



(2) 



uj - £ - Un- a — Y, a (uj) + ir 

where n a denotes the electron number in the dot and 
T, a (oj) is the self-energy to be calculated in terms of the 
zeroth-order Green's function G^\u) — (u> — Eq + ir)~ 1 , 
r/2 = T L = T R = irp c (0)V 2 . p c (0) denotes the den- 
sity of states (DOS) of conduction electrons at the Fermi 
level. In our modified perturbation theory, we intro- 
duce an effective energy level io a in place of the Hartree- 
Fock level £o + Un- a and determine £ UC r to satisfy the 
Friedel sum rule exactly. We term this method the self- 
consistent energy, level second-order perturbation theory 
(SCEL-SOPT)lj Thus, the effective first-order Green's 
function is given by 

GW( W ) = J——. ( 3 ) 



u> ■ 



ir 



Using this Green's function, we first calculate the or- 

(2) 

dinary second-order self-energy (w). Then, we in- 
troduce the modified self-energy which is correct ixL the 
atomic limit ( T — > , p a {ui) — > 5(uj - £o CT ) ) asEI'Ei 



where 



{7(1 - n { l' a ) + £ - soo 



B = 



n[c; 



and 
n[G a ] = 



dw/(a;)(-l/7r)ImG (T (w) 



(4) 



(5) 



(6) 



Note that i? vanishes in the electron-hole symmetric case 
( e =-U/2 ). 

Next, we construct the second-order Green's function 

as 

1 



Eq — Un 



(i) 



s«r(w) + ir 



(7) 



Using this G^ ; (w), we calculate the second-order elec- 
tron number n {2) = n[G { a ] ]. Further more, Un_ a in the 

(2) 

denominator of eq. (7) is replaced by Un_ a to obtain 
the solutions for £o over a wide range. We calculate 
— n[G^] again and ?ii 2 ' )FS = n FS [Gi 2 ' ) ] from the 
Friedel sum formula 



(8) 



Then, we determine £Q a so as to satisfy the relation 
— 71ct 2 ' F ^ (Friedel sum rule). 
The present scheme (SCEL-SOPT) is applicable only 
to zero temperature, since the Friedel sum formula is 
only valid in that case. However, we extend this scheme 
approximately to investigate the properties of a dot at 
finite temperatures by using Eq^ determined at zero tem- 
perature.E3 

To confirm the accuracy of our approximation, we 
compare several properties with those of the Bethe Ansatz 
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Fig. 1. The spin susceptibility multiplied by the temperature is 
shown for U=4 , T=l and eo=— 2. The full line represents the 
results calculated by the Bethe Ansatz solution and the circles 
represent our calculated results. 
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Fig. 2. The specific heat is shown for U=4 , T=l , h=0 and 
£0=— 2. The full line represents the results calculated by the 
Bethe Ansatz solution and the circles represent our calculated 
results. 



solution to the impurity Anderson model. Figure 1 shows 
Txm where \m is the spin susceptibility for U=4, T=l 
and £q=— 2. The unit of energy here and throughout 
the present letter should be meV for typical quantum 
dots. Th£ full line represents the results of the exact 
solutionE-$ and the circles represent those of our calcula- 
tions. At high temperatures, the exact solution gradually 
approaches 1/8. Our result shows qualitative agreement 
with the exact solution over a wide range of tempera- 
tures. Figure 2 indicates the specific heat C — dEtot/dT, 
where i?tot is the total energy calculated from 



ot 



dw/(«)(-i) 



xlm[{c- -{Un™+^{w))}GM{w)). (9) 

Here, Ecr(w) is the modified self-energy given by eq. (4). 
The global feature of C(T) is well reproduced by the 
present calculation. 

As we have assumed a single energy level £o in a quan- 
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Fig. 3. Temperature dependence of conductance is shown for 
[7=4 , h=0 and r=l(a), 0.5(b). 



turn dot, the current takes the for nil) 

1 = \ E / - M<*)}r'p*W, (io) 

where r' = r L r R /(r L + T R ) and f v (u) is the Fermi 
function of the lead v. We utilize eq. (10) to obtain 
the differential conductance (dl/dV) at zero bias voltage 
which is equivalent to the linear conductance, where the 
bias voltage is given by the difference of the chemical 
potential between the left and right leads as V=hl—Pm- 
We can use an equilibrium DOS for p a (^) in eq. (fO) 
to calculate the transport properties in a linear response 
regime. 

Figure 3 shows the temperature dependence of the lin- 
ear conductance through a quantum dot for U—A, h=0 
and r=l(a), 0.5(b) as a function of the gate voltage, 
where h = ^gfiBH. The full, broken, dotted and dash- 
dotted lines represent the results for T=0, 0.05, O.f and 
0.2, respectively. At low temperatures, the Kondo peak 
of DOS emerges at the Fermi level in the symmetric case, 
and hence, the conductance takes the value 2e 2 /h due to 
the Kondo-assisted tunneling in that case. In the case 
of electron-hole asymmetry, the Kondo peak shifts from 
the Fermi level and the conductance is reduced. As the 
Kondo effect does not appear in the extreme asymmet- 
ric case, So <— 4 or £o >0 for example, the conductance 
does not depend on the temperature. In Fig. 3(b), the 
conductance exhibits two peaks for T=0.2 which is noth- 
ing but the Coulomb oscillation. The peaks appear when 
the discrete energy level of a quantum dot Eq or Sq + U 
agrees with the Fermi level. These two peaks are not 
seen clearly in Fig. 3(a) since the width of the reso- 




Fig. 4. Magnetic field dependence of conductance is shown for 
(7=4 , T=l and T=0(a), 0.05(b) and 0.1(c). 



nant level T is large. These temperature dependencies 
of the conductanrj^were already obtained theoretically 
by other groups,BE3ES ) and our numerical calculation is 
consistent with their results. Recent experimentsEm' n 
have revealed Kondo-assisted tunneling at low tempera- 
tures that demonstrates qualitative agreement with our 
results. 

Figure 4 shows the magnetic field dependence of the 
conductance for U—4 and T=l. The full, broken and 
dash-dotted lines represent the solution for h=0, 0.2 and 
0.3, respectively. We have included an external mag- 
netic field H by shifting the dot level as £o ± h for up 
and down spins. The chemical potential /^„'s are as- 
sumed to be unchanged by the field. Under the magnetic 
field, the Kondo peak shifts to the right and left for up 
and down spins, respectively, by Zeeman energy. Then 
the DOS at the Fermi level, which determines the trans- 
port properties, decreases, and hence, the conductance 
near electron-hole symmetry is reduced. This behavior 
appears to be-consistent with that recently observed in 
experimens.LTtT The above feature is similar to the case 
of reduction in conductance with increasing temperature 
in Fig. 3. In the extreme asymmetric case, the oppo- 
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site feature emerges which enhances the conductance by 
the application of the magnetic field, namely, positive 
magnetoconductance. As regards the sign of the magne- 
toconductance, we assume that, in the range of negative 
magnetoconductance, the dot contains an odd number of 
electrons and the Kondo-like effect occurs. Otherwise, 
there are an even number of electrons and the Kondo 
effect does not occur. Furthermore, we observe that in- 
creasing the temperature weakens the magnetic field ef- 
fect. 

In conclusion, we have studied the transport properties 
through a quantum dot in the Kondo regime by the mod- 
ified perturbation theory (SCEL-SOPT), which modifies 
the self-energy so as to reproduce the correct atomic limit 
and to fulfill the Friedel sum rule exactly. We extended 
this theory to study the finite temperature properties of a 
quantum dot. The accuracy of our scheme was confirmed 
by comparing the spin susceptibility and the specific heat 
with the exact solution. It has been demonstrated that 
our scheme can be used for a wide range of tempera- 
tures and magnetic field strengths. As regards the trans- 
port properties through a dot, we noted the reduction 
of the conductance near electron-hole symmetry under 
magnetic field at low temperatures, which explains the 
recent experimental resultsiM We also showed the posi- 
tive magnetoconductance in the case of asymmetry. Such 
properties may also be observed in the experimental re- 
sults. During the preparation of the present manuscript, 
we discovered a paperEf which also calculated the mag- 
netic field dependence of the conductance through a dot 
by QMC. The result is consistent with that of our calcu- 
lation near the electron-hole symmetric case. However, 
because of the use of QMC, their calculation does not 
cover sufficiently low temperatures or a wide range of 
gate voltages. Our approach is much simpler with suffi- 
cient accuracy and may have wide potential applicability 
to various situations in quantum dot systems. 

The computation was performed using FACOM VPP500 
in the Supercomputer Center, Institute for Solid State 
Physics, University of Tokyo. This work is supported by 
a Grant-in- Aid for Scientific Research No. 11640367 from 
the Ministry of Education, Science, Sports and Culture. 
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